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Abstract 

Ultrafast photoinduced dynamics of electronic excitation in molecular dimers is drastically af- 
fected by the dynamic reorganization of inter- and intra- molecular nuclear configuration modeled 
by a quantized nuclear degree of freedom [Cina et. al, J. Chem Phys. 118, 46 (2003)]. The dy- 
namics of the electronic population and nuclear coherence is analyzed by solving the chain of cou- 
pled differential equations for population inversion, electron-vibrational correlation, etc. [Prezhdo, 
Pereverzev, J. Chem. Phys. 113 6557 (2000)]. Intriguing results are obtained in the approximation 
of a small change of the nuclear equilibrium upon photoexcitation. In the limiting case of resonance 
between the electronic energy gap and the frequency of the nuclear mode these results are justified 
by comparison to the exactly solvable Jaynes-Cummings model. It is found that the photoinduced 
processes in the model dimer are arranged according to their time scales: (i) fast scale of nuclear 
motion, (ii) intermediate scale of dynamical redistribution of electronic population between excited 
states as well as growth and dynamics of electron-nuclear correlation, (iii) slow scale of electronic 
population approach to the quasi-equilibrium distribution, decay of electron-nuclear correlation, 
and decrease of the amplitude of mean coordinate oscillation. The latter processes are accompa- 
nied by a noticeable growth of the nuclear coordinate dispersion associated with the overall nuclear 
wavepacket width. The demonstrated quantum relaxation features of the photoinduced vibronic 
dynamics in molecular dimers are obtained by a simple method, applicable to systems with many 
degrees of freedom. 

PACS numbers: 02.30.Jr, 05.10.Gg, 31.50.Gh, 82.20.Kh, 89.30.Cc 
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I. INTRODUCTION 



Recent achievements of femtosecond spectroscopy^ open new windows into the world of 
electron and exciton transfer dynamics in molecular systems, which play dominant roles in a 
variety of problems in physics, technology, atmospheric photochemistry, and biology^'^'^'^'^. 
Electron and exciton transfer drive formation and breaking of chemical bonds, corrosion 
reactions, ion tunneling microscopy, enzymatic activity in living cells, energy storage by 
adenosine triphosphate and harvesting of solar energy. For instance, the harvesting and 
storage of solar energy in the natural photosynthetic complexes, such as bacteriochlorophyll^ 
or Rhodobacter sphaeroids^, rely on the specific alignments of the light-harvesting units that 
favor efficient light absorption and subsequent energy transfer. The natural and artificial 
light harvesting complexes (LHC)® involve exciton transfer of the primary photo-excitation 
in the antennae complex^, followed by a subsequent electron transfer that takes place in the 
reaction center of the natural LHC and porphyrin arrays of artificial LHC^°'^^. 

The dynamics of electron and energy transfer, and formation of entangled electron-phonon 
states are investigated by a wide arsenal of femtosecond experimental techniques, which 
by application of one or several short laser pulses, provide time resolved information on 
couplings, correlations and states of molecular aggregates^^'^^'^^'^^'^^'^^. The multi-pulse time 
resolved techniques cover various frequency ranges, including X-ray, visible and infrared, 
and serve for determination of population, location and phase of electron-nuclear states of 
molecular systems. Thus, a recently developed Raman X-ray spectroscopy follows electronic 
state dynamics^^. Laser pulses in the infrared diapason are used to investigate populations 
and mutual correlations of vibrational modes^^'^*^. Visible frequency short laser pulses are 
applied to measure electronic state populations and state-to-state correlations^^. Positions of 
nuclear wavepackets are detected by pump-probe spectroscopic techniques^ . Relative phases 
of electron-nuclear wavepackets of molecular aggregates are obtained using the nonlinear 
wavepacket interf eromet ry , 23 , 24, 25 _ 

Theoretical methods describing exciton and electron transfer dynamics appeal, in one way 
or another, to both electronic and nuclear degrees of freedom^^'^^'^^ and most commonly in 
chemical physics employ the notion of a reaction coordinated^. The reaction coordinate 
quantifies changes in the electronic states due to nuclear mode dynamics and is based on 
the fact that transfer is most effective for specific configurations of the nuclear subsystem. 
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These nuclear configurations correspond to crossings of tlie potential surfaces and match 
the Franck-Condon window for the source and target states. In many cases, particularly for 
large systems, the nuclear subsystem leaves and never returns to the optimal configuration 
region due to stochastic noise and destructive interference among nuclear modes of different 
frequencies involved in the dynamics. This dynamic mismatch and dephasing of the nuclear 
modes is responsible for relaxation and decoherence^" in the exciton and electron transfer 
dynamics. 

The coupled electron-phonon dynamics can be described by a number of theoreti- 
cal approaches developed in chemical physics^^, including the Gaussian wave-packet ap- 
proaches^^'^'^, the semiclassical approximations^^ based on Feynman path integrals^^, and 
non-adiabatic molecular dynamics^^''^^'^^. The quantum-classical mean- field approach and 
its multiconfiguration mean-field and surface hopping extensions are applied to a wide range 
of problems, including gas-phase scattering, surface^, solution and biological chemistry^^. 
Semiempirical approaches introduce quantum corrections to classical mechanics. The re- 
cently developed Quantized Hamilton Dynamics (QHD) method^°'^°'^^'^^'^^'^^ offers a reduced 
description of quantum nuclear motion by disentangling coordinate, momentum, dispersion 
and higher order variables for each nuclear mode. A many body closure"^^ terminates the 
infinite chain of coupled equations reducing the number of equations and saving the calcu- 
lation time. While the semiclassical and mean-field approaches are very effective from the 
numerical point of view, they provide approximate solutions. In particular, the mean-field 
approaches typically cannot describe the dynamic creation of superposition and entangle- 
ment between states leading to branching of nuclear dynamics. 

Theoretical models that are able to describe the superposition and branching feature of 
coupled electron-phonon dynamics are available in the solid-state theory^^ and quantum op- 
tics theories of atom-field interactions^^. The quantum optics approaches to the description 
of photon modes coupled to atomic quantum states^^ can be adopted to study electron- 
phonon states in chemical, biological and nano-systems. A simple eigenstate solution for the 
correlated dynamics of two electronic states of an atom (fermion degree of freedom) coupled 
to a laser mode (boson degree of freedom), presented in 1963 by Jaynes and Cummings^^'^°, 
takes an important place among the optical methods. The Jaynes-Cummings solution was 
further developed by the operator algebra methods^^. Known are the operator solutions 
for the lowering and raising operators^^ and for the electronic state population^^ that shows 
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collapses^^ and revivals^^. Interesting expressions are found for the uncertainty relationships 
for the operators describing the two- level atom interacting with a photon field^^. The gener- 
alized operator solutions lead to a well-developed formalism in the Fokker-Planck equation 
formulation^^. The collapses and revivals in the two-level atom excitation dynamics of the 
Jaynes-Cummings model represent the simplest case of quantum relaxation common for 
exciton and electron transfer in molecular complexes^^'^^'^°. 

The quantum optics and QHD approaches are combined in this paper for the description 
of molecular aggregates, whose photoexcitation dynamics is strongly affected by dynamic 
rearrangements of vibrational degrees of freedom. The combination of methods is very 
effective for the calculation of the relaxation behavior of exciton transfer in molecular ag- 
gregates. Extending the original QHD approach that focuses onto semiclassical dynamics 
of nuclear modes, a QHD approximation is developed in present for the dynamics of the 
coupled electronic and nuclear variables, concentrating on the difference of the electronic 
state populations, i.e., population inversion. The resulting equations are very simple to be 
applied to large condensed phase chemical systems. Without reference to a thermal bath, 
as typical in the quantum relaxation theory, the QHD approximation for the dynamics of 
the Jaynes-Cummings model shows quantum relaxation features. The relaxation charac- 
ter of the exciton and electron transfer dynamics arises naturally in our approach due to 
destructive interference of quantum states, technically similar to quantum beats. The com- 
plex dynamics of the electron population in the Jaynes-Cummings model is made more 
understandable by a time-scale hierarchy^^'^^ of the relevant dynamical processes. 

The paper is organized as follows. The relevant theoretical tools are introduced in Sec- 
tion II. The calculated quantum dynamics are analyzed in Section HI, including a detailed 
discussion of the vibronic wavepacket. Section IV considers approach of the system to a 
state with large uncertainty. The physical processes in the system are organized in accor- 
dance with their characteristic times in Section V. The scenarios of experiments that can 
be described by the current method, and where the dynamical features found in this paper 
can play a dominant role, are discussed in Section VI. Finally, Section VII summarizes the 
main conclusions of this work. 
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II. MODEL 



Consider dynamics of electronic states in a molecular system. Typically, the initial pop- 
ulations of the states are derived from a thermal equilibrium. An external perturbation 
such as a femtosecond optical pulse quickly, within 10^^^ — 10^^^ s breaks the equilibrium 
by enhancing population of excited states. The non-equilibrium populations evolve in time. 
The populations of the excited states are redistributed within an intermediate time interval 
of a few picoseconds, 10~^^s. The population of the electronic states dephases due to in- 
teraction with inter and intramolecular vibrational modes. On a relatively long time scale 
of nanoseconds, 10~^s, the molecular system is de-excited by spontaneous emission induced 
by interaction with zeroth order vacuum oscillations of photon modes, or by some other 
mechanism. Upon electronic de-excitation, additional electron-vibration dynamics returns 
the system to thermal equilibrium. 

The time-dependent processes that occur shortly after pumping are of interest in the 
present work. The focus is on dynamics of a pair of single-excitation localized electronic 
states z = 0, 1. One of the states z = 1 is pumped initially and is coupled to and exchanges 
population with the second state. The two electronic states will be described in terms of the 
second quantization operators, creating cf , or annihilating q the excitation of localized state 
i. The Hamiltonian describing the picosecond t ~ lO^^^s dynamics of the excited states can 
be written in terms of the creation and annihilation operators as 

Ho = ^eicfci + ^^^Jijcfcj (1) 

i i j 

The first term contains the electronic state energy times the operator of number of quanta, 
whose expectation value equals one (c^q) = 1 if state i is populated and zero if state i is 
empty. The second term describes the coupling between the electronic states by removing 
population in state j, Cj and creating population in state i, cf with the coupling constant 
Jij. In general, both energy and coupling depend on nuclear configuration. For small 
fluctuations of the nuclear coordinate q the dependence of energy and coupling on q can be 
well represented by the zeroth and first order terms in the Taylor series expansion 

e^{q) = + z = 0, 1. (2) 

Generally, the expansion coefficients are specific for each electronic state. The Taylor ex- 
pansions of the coupling constants may also be different for the forward and backward 
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transitions 



7 _ t(0) ,^9Joi I 



^10 = jl?+g^U + .... (3) 

Typically, however, the coupling constants are identical for both directions Joi = ^lo, 
dJoi/dq = dJio/dq. The position independent terms in the expansion of the coupling 
constant can be eliminated by diagonalizing of the g-independent part of the electronic 
Hamiltonian. 

The evolution of the nuclear coordinate q augments the original electronic Hamiltonian (1) 
with a vibrational term, which equals /{2m) + muo'^q^ /2 in the harmonic approximation. 
The vibrational Hamiltonian is the same for both electronic states. Multiplication of the vi- 
brational Hamiltonian {p^ / (2m) + muj'^q^ /2) Id by a unit operator in the electronic subspace 
lei = X]j ^t^i leads to the following electron-phonon Hamiltonian 

- - E ^t. { ^ . , . 14 . ? ^ ^- {- 1 } 

Consider the coordinate-dependent terms indicated in (4) by the underbrace sign in more 
detail. The equilibrium position for the harmonic motion correlated with the electronic 
state i can be obtained explicitly by completing the square: 

/ muj'^ 9 muj'^ \ I 2 , f ^'i \^ / e!- ^ ^ 
2 2 muj'^ V mcu^ / \ mcu^ 




where = dei/dq, Ei = ef'' + (e^)^/2mu;^, is the renormalized state-specific electronic 
energy and = e[{muj'^)"^ . The electron-phonon Hamiltonian takes the form 



H = ^c-c, u ,.A + y: y: etc, L . II 



(6) 



Without loss of generality go = 0, gi 7^ 0. The Hamiltonian describing two electronic 
states coupled through a vibrational mode contains two harmonic potential energy surfaces 
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corresponding to the two electronic states and centered at and qi. The states are coupled 
via a coordinate-dependent term similar to the Holstein Hamiltonian^^ 



H 



Co Co 



_ lower potential 



El + rauP' jl (g — gi)^ /2m 



upper potential 



r + +1 



(7) 



transfer 



A two- level system of single excitation states is conveniently represented as a spin-1/2 par- 
ticle in terms of three spin-projection operators ^x, Sy, S^, which can be expressed us- 
ing the lowering and raising operators = 1/2 (5+ + S-), Sy = i/2 {S+ — S-), = 
1/2 {S+S--S-S+). 



S±S^ = 1/2 ±5„ 



(8) 



The unique relationships between the spin-1/2 raising and lowering operators and the cre- 
ation and annihilation operators of the electronic system in the language of second quanti- 
zation are given by 



Co Co = 


|0) (0| 




c|ci = 




= S-^S ^ 


cfco = 


|1)(0| 


= s^^ 


CqCi = 


|o)(i| 


= s_. 



(9) 

The evolution of the vibrational coordinate and momentum relevant for experiments 
occurs on an atomic scale demanding a quantum mechanical treatment. The quantum 
coordinate q and momentum p operators are conveniently replaced by the harmonic raising 
a"*" and lowering a operators 



p 



A/ — ^ (a+ + a) 
V 2muj ^ ' 

hmuj , 

V 2 ^ 



— a) 



(10) 



The above operator transformations translate the problem into the second quantization 
language of raising and lowering operators for both the electronic and vibrational degrees 



of freedom. The relevant dynamics will be represented by product operators containing 
spin and harmonic raising and lowering terms. The product operators describe resonance 
exchange of quanta between vibrational and electronic subsystems. The shift in the equilib- 
rium coordinates of the two harmonic potentials is represented in the second quantization 
language by the dimensionless parameter a = q^ymuj/2h + ip^/h/2muj, related to the orig- 
inal coordinate shift 

The zero of electronic energy is chosen so that the initially unoccupied lower state i = has 
zero energy. The physically relevant electronic energy gap HQ = Ei — Eq is independent of 
energy origin. With all above changes, the electron-phonon Hamiltonian reads 

H = hujS^S+ [a+a + 1/2] + huoS+S. [{a+ + a*){a + a) + 1/2] 

+g {S+ + 5_) (a+ + a) + S+S^ [hn - huja*a] (12) 

where g = ^Jfijlmuj dJ /dq. 

The Hamiltonian (12) represents a dimer with a pair of electronic states |0) and |1). The 
states are separated by energy difference Ml and are coupled to a single harmonic vibrational 
mode. The equilibrium positions of the potential energy surfaces describing the vibrational 
mode in the two electronic states are shifted with respect to each other by gi, specified by 
the dimensionless nuclear reorganization parameter a, as illustrated in Fig. IThis model 
is known as the Marcus model^^ that has an enormous range of applications to exciton, 
electron, proton transfer and many other chemical reactions^^'^^''^^. 

In the limit of small nuclear reorganization the electron-phonon Hamiltonian transforms 
into a simpler form, known in quantum optics as the Jaynes-Cummings Hamiltonian^^'^*^. 

i/jCM = ^ {a+a + 1/2) + nS+S^ + g {a+ + a) {S+ + S.). (13) 

Here and below h = 1. The Hamiltonian contains three terms corresponding to the electronic 
subsystem, vibrational subsystem, and excitation transfer. The approximations below are 
developed for the Jaynes-Cummings Hamiltonian, but can be equally obtained for a non- 
adiabatic Hamiltonian with the momentum dependent coupling g (a~^ — a) (5*+ + S'_) rather 
than the position dependent coupling g {a+ + a){S+ + S^). A further simplification of the 
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Jaynes-Cummings Hamiltonian is obtained by invoking the rotating wave or resonance ap- 
proximation that keeps only those couphng terms that preserve the total number of quanta 
in the combined electron-vibrational system. 

HjcM = ^ [a+a + 1/2) + nS+S^ + g (a+5„ + aS^) . (14) 

The dynamics that follow from this final form of the Hamiltonian are investigated below. 

The expectation value of the operator is the main focus of the present study. This 
expectation value gives populations of the electronic states of the dimer and, more generally, 
is related to the value of a reaction coordinate in exciton, electron, proton transfer and other 
chemical processes. In the Heisenberg representation of quantum mechanics, the expectation 
value of Sz evolves in time through the time-dependence of the operator. The wave function 
remains fixed and is specified by the initial conditions. The Heisenberg equation of motion 
(EOM) for the operator Sz is given by its commutator with the Hamiltonian (14). The 
time-derivative of depends on other time-dependent operators, whose EOMs are also 
obtained by commutation with (14). This leads to an infinite hierarchy of EOMs. Noting 
that Sz = 1/2 {S^S^ — S^S^), the infinite hierarchy can be written in a compact form using 
only three additional operators 

a = a~^S- + aS^, 
P = a~^S^ — aS+, 

7 = a+a + S+S. + (15) 

that involve the non-interacting Hamiltonian 7, the interaction part of the Hamiltonian 
a, and the auxiliary operator /3. The following hierarchy of Heisenberg EOMs is obtained 
starting with the expectation value of the population inversion operator Sz 

'Jt ^"^"^ = ^'^^"^ ' 

^|(/57") = -5(«7")+£(£^y;^, 

needs closure 

^J^ {S.ll = 9 {Pll (16) 

with n = 0, 1, 2, . . .. The detuning 6 = Q — u denotes the difference between the electronic 
energy gap Q and the vibrational frequency u. The Heisenberg EOMs for the operators are 
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arranged in blocks of n-th order with three equations per block. The evolution of the lower 
order block depends on the evolution of the higher order block through the coupling term 
labeled in (16) by underbrace. 

Useful approximations to the exact solution of the infinite chain of equations (16) can be 
obtained by limiting the number of equations to a few lower order blocks and decomposing 
the higher order coupling term into a product of lower order terms in the spirit of the 
quantized Hamiltonian dynamics (QHD) approach originally introduced for the classical- 
like expectation values of the position and momentum operators^'''^'^. Decomposition of 
the higher order expectation values into products of the lower order ones terminates the 
infinite hierarchy (16) within a finite number of equations. With no a priori knowledge of 
the relationship between the expectation values of the higher and lower order operators, 
the contributions of all possible products of the lower order operators that represent the 
higher order operator are taken with equal weights^^'^'^. The closure applied to decompose 
higher order moments of the position and momentum operators into products of the first 
order expectation values gives classical Hamiltonian mechanics. Decomposition of the higher 
order moments of position and momentum into products of the first and second moments 
leads to the well known Gaussian approximation^^. Generally, the QHD procedure allows 
one to obtain simple classical-like EOMs for essentially quantum mechanical characteristics, 
such as zero-point energy, tunneling, quantum state populations, state-to-state correlations. 

In present, the hierarchy (16) is terminated at n = 1 requiring decomposition of the Szj'^ 
term. Using the general expression for the decomposition of a triple product^° 




(17) 



this term is decomposed into 



(S.l) (7) + (S.) {l') - 2 (S.) (7)^ 



(18) 
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yielding a closed system of six differential equations 
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Usually, closures produce systems of non-linear differential equations, such as the non-linear 
classical Hamiltonian EOMs obtained from the Heisenberg quantum EOMs by decomposi- 
tion of the expectation values of higher order moments of position and momentum operators 
into products of the first order moments^'''^''. The special properties of the Jaynes-Cummings 
Hamiltonian in the rotating wave approximation, [7, H] = in particular, make Eqs. (19) 
into a system of linear differential equations. The operators 7 and 7^ appearing in the clo- 
sure (18) are integrals of motion of the Hamiltonian (14). The quadratic (7)^ term in (19) is a 
constant specified by the initial conditions. Equations (19) give a very simple approximation 
to the dynamics of the electron-phonon system. 



III. ANALYSIS OF DYNAMICS 



The closed system of Heisenberg equations (19) can be solved analytically for the ex- 
pectation value of the population inversion {3^). The time-evolution of {3^) is given by a 
superposition of two cosine functions, 

{3^} (t) = {3^} (t = 0) H ^ (1 - coscjit) 

Tr2 (1 - cos UJ2t) , (20) 

ZUJ2 

whose frequencies are defined by 

ujI = 5^ + 2g^ (70 + y^) , 



2 



6' + 2g' (70 - ^0) 



70 = (7) {t = 0). (21) 

The difference in the squares of the frequencies is given by the product of the square of the 
electron-phonon coupling constant and the square root of the mean number of quanta 
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in the system ^/tq. As should be expected, a perfect oscillation in the electronic state 
population is observed with zero electro n-phonon coupling. As the strength of the electron- 
phonon coupling grows, the difference between the two frequencies increases, leading to the 
dephasing of the electronic state population and a corresponding decrease in the oscillation 
amplitude. Fig. 2. The oscillation of population inversion will proceed more rapidly if the 
vibrational mode contains more energy. Vibrational motion of a large amplitude accelerates 
transfer of electronic population. The time evolution of the coupled electron-phonon system 
is illustrated in Fig. 2for the initial conditions 

{l)\t=0 = < >' 

(«)lt=o = 0, (/3)U = 0. (22) 

The vibrational mode is prepared in a quasiclassical coherent state described in the coordi- 
nate representation by a Gaussian displaced from the vibrational equilibrium. All electronic 
population is localized in one state. There is no correlation between the electronic and vi- 
brational subsystems at the initial time. The evolution predicted by Eq. (20) is compared 
with the exact solution^^. The exact and approximate solutions coincide with good precision 
until 27rgt < 1 holds^^, corresponding to t < 6.366 and t < 0.6366 in Figs. 2(a) and 2(b), 
respectively. The approximate solution is analytic, in a simple closed form, Eq. 20 compared 
to the exact solution involving an infinite series summation^^. Note that within the time 
interval determined by 27rgt < 1, the dynamics of inversion has relaxational character and 
usually fits the expression ^ cos (2y/^gtj exp {—g'^t^Y^. 

The sum of the cosine functions in Eq. (20) forms beats. The fast oscillation of the 
population inversion in Fig. 2 gives excitation transfer between the electronic states. The 
envelope of the beats describes relaxation in the oscillation of the electronic population 
transfer due to the vibration induced dephasing. The inversion and relaxation tZ^ rates 

^+ " 2 ' 

= (23) 

depend on the number of phonons. The sum and difference of the frequencies, determine the 
beats of the electronic population in Eq. (20). The analytic dependence of the rates of in- 
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version (faster component) and relaxation (slower component) on the vibrational amplitude 
(70) and electron-vibrational detuning 5 = u; — ^2 is illustrated in Fig. 3. 

The dependence of t^^ on S and (70) becomes clear from the following series expansion: 
For (5^1 and (70) ^ 1 the frequencies t±. of the population inversion oscillation and 
relaxation can be expanded in powers of 6 and 1 / a/ (70) 

The inversion rate t^'^ has minimum at 5 = 0, the relaxation rate tl^ has maximum at 5 = 0. 
A larger detuning leads to faster oscillations of population and a slower relaxation. 

The dependence of t± on the amount of vibrational energy is considered in the 5 = 
limit, which has been explored in Fig. 2. Expansion of the sum and difference frequencies 



up to the third order in 1/ a/ (70) with 6 = takes the form: 
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- -7= f l + ^^l +•••• (25) 

The frequency of population inversion grows with system energy (70) = a*a, already in 
the zeroth order of expansion. The relaxation frequency decreases as inverse of system 
energy, Fig. 3. Considering the vibrational mode as a quasi heat bath with respect to the 
two-level system, it is intuitively expected that hoter bath yields quicker relaxation. The 
decrease of the relaxation rate with increasing phonon energy occurs since the oscillator is 
not in a thermal state, but approaches at large (70) the classical regime. A classical oscillator 
coupled to a two state electronic system yields oscillation in the electronic population without 
relaxation. For small (70) < 1 the oscillation and relaxation frequences coincide. For large 
(70) > 1, as indicated by expansion (25) the oscillation frequency grows while the relaxation 
frequency decreases. The frequencies display branching with energy growth. Fig. 3. 

The approximate solution Eq. (20) not only gives the first few oscillations of the electronic 
population corresponding to several forward and backward population transfer reactions, but 
also well reproduces the overall dephasing envelope that is associated in the thermodynamic 
limit of many vibrational modes with relaxation to equilibrium. It is quite remarkable that 
this quasi-equilibration is observed transiently already with a single vibrational mode! 
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It is instructive to look further into the dynamics of the population inversion (5'z). The 
dynamics of this observable is conjugate to the dynamics of the expectation value of the 
correlation operator defined in Eqs. (15). As shown by the exact solution in Fig. 4, 
after a few oscillations the population inversion stalls at zero indicating that the electronic 
population is equally delocalized between the two states. This situation is referred to as 
"collapse" After a while, the oscillations reappear, Fig. 4, producing a "revival" It 
is important to understand in more detail this type of dynamics, in particular, the fate of 
the energy and phase dynamics during the "silent period", after the collapse and before 
the revival. The electronic contribution to the total energy is stored in an unusual form 
during the silent period. The electronic population is distributed equally between the two 
states. This situation is equivalent to the high temperature limit of the populations of a 
two- level system coupled to a thermostat. Here the nontrivial dynamics appear due to the 
entanglement of the electronic and vibrational degrees of freedom. 

It is intriguing what happens to the phase of the oscillation during the silent period. The 
phase is preserved in a hidden manner to reappear later when the silent period is over. Fig. 4. 
It may be expected that a simple operator or a combination of several operators exist, whose 
expectation value maintains the oscillation, while the oscillation of the population inversion 
stalls. The dispersion 

aA = {{A-{A)f). (26) 

of the vibrational coordinate {A = q) is the desired expectation value. The vibrational 
dispersion starts at the minimal value at time zero, grows with time and maintains the 
oscillation during the silent period. Fig. 5. 

The origin of the oscillation in the dispersion of the vibrational coordinate is illustrated 
by the evolution of the vibrational wavepacket in Fig. 6 that shows wavepacket snapshots 
at several times. Initially, the wavepacket is Gaussian. As time evolves, the electron- 
vibrational coupling splits the vibrational wavepacket into two branches correlated with the 
two electronic states. The dispersion of the vibrational coordinate represents the width 
of the overall wavepacket. The width oscillates as the two wavepacket branches evolve 
nearly independently as determined by their own electronic states. The width of the overall 
wavepacket is determined by the separation between the branches. The width is maximal 
when the two wavepackets are far apart and minimal when the two wavepackets penetrate 
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through each other. After the silent period when the population inversion resumes its 
oscillation, the wavepacket width oscillation decreases in amplitude, Fig. 6. 



IV. ANALYSIS OF UNCERTAINTY 

Dephasing of the electronic subsystem due to coupling to vibrations is ubiquitous in 
photochemistry. The silent period when the dynamics of the electronic population in the 
Jaynes-Cummings model temporarily stops due to the vibrationally induced dephasing is 
investigated further. Consider the uncertainties of the observables describing the vibrational 
and electronic modes. The vibrational wavepacket with two peaks that exist during the 
silent period. Fig. 6, has a large uncertainty. As shown before, the phase of the electronic 
subsystem is preserved in a hidden form by the uncertainty of the vibrational coordinate. 

The uncertainty relationship for dispersions cr^, a^, Eq. (26), of Hermitian operators A 
and B is commonly written in a product form as^^ 



■ 0"B > 



(27) 



An alternative relationship between the dispersions in a sum form can be obtained based 
on the positivity of the norm of operator F = A — (A) -\- i (B — (B)), {F*F) > — >• 
{{A - {A)f + {B- {B)f + i[A-{A),B- {B)]) > 0, and 

aA + (Tb + I {[A, B]) > 0. (28) 

The sum + as + i {[^,B]) gives a good second-order combination of operators, whose 
expectation value can be included into the chain of Heisenberg EOMs. In application to 
the vibrational subsystem, the sum uncertainty relationship (28) means that the sum of 
dispersions for the vibrational position and momentum 

A = + a = q\/ 2muj / h 



must be larger than zero 



B = i{a^ -a)= p^2/{hmuo) (29) 



(a+a) - (a+) (a) > (30) 
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Consider the dispersion of coordinate in more detail. It can be written in terms of the 
raising and lowering operators in the form 

= (31) 

= (\a~^ + a- (a+ + a)] ' 

h 



{ (aa+) - (g) (a+) + (a+a) - (a+) {a) + (a+)'^ - (a+)' + (a^) - {af } . 



2mu ■ 

^ ^ ' 

cross term i jj i 

ladder term 

The coordinate dispersion includes two kinds of terms, the ladder dispersion terms involving 
operator squares that increase or decrease the number of quanta by two, and the cross terms 
including products of conjugate ladder operators that preserve the number of quanta. The 
mean value of the ladder term oscillates quickly and gives fringes in Fig. 5. The cross term 
maintains the rotating wave approximation preserving the number of quanta, and evolves 
slowly. Similar decomposition applies to the dispersion of momentum 



(Jr, 



2 



{[a+ - a - {a+ - a)] ) (32) 
(aa+) - (a) (a+) + (a+a) - (a+) (a) - [((a^)') - (a^)' + (a'> - (a)' 



cross term i jj ^ 

ladder term 



The difference from the position operator dispersion is only in the negative sign by the 
ladder term. The quickly oscillating ladder terms cancel out in the sum of the coordinate 
and momentum dispersions, leaving only the cross terms 

i + = i^-. + = I (fM + V2 - (oW } = . (33) 

V. total energy semiclassical / 

The slowly oscillating sum of dispersions (33) is nothing but a purely quantum contribution 
to the vibrational energy defined as the difference between the total vibrational energy 

i^t^otai = ^(/) + = Ma^^ + 1/2), (34) 

and the semiclassical contribution to the total energy 

^stmi-ci = ^{P? + "^{q? = nuj{a+){a), (35) 

formed of the first order classical-like averages (g) and (p). The purely quantum contribution 
to the vibrational energy provides the envelope of the coordinate dispersion in Fig. 5. Sub- 



stitution of the commutation relation [p,q] = ih [-^^B] = q^j2muj /h,p\j2/muh 



2i 
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and the sum of dispersions (33) into the general formula (28) gives the uncertainty relation- 
ship (30) for the vibrational dispersions. 

Next, consider the sum uncertainty relationship, Eq. (28), for the electronic, generally 
fermionic subsystem with 

B ='-{S^- 5_) = S,. (36) 



The inequality (28) gives 



- (5+) > 0. (37) 



Similarly to the vibrational subsystem, the dispersions of the electronic S'x and Sy operators 
can be decomposed into the cross and ladder terms. The dispersion (Ts„ reads 

as^ = l([S+ + S„-(5+ + 5_)f> (38) 

= - { {S+S^) - (5+) + - (5+) + (Si) - {S+f + {S'_) - {S^f } . 
4 V — . ^ — ^ — . ' 

cross term ladder term 

Since the second order averages in the last equation add to a constant, {S^S^ + S^S+) = 1 
by completeness, and (5*+)^ = = due to finite dimensionality of the electronic basis, 

the electronic dispersion involves only the first order averages 

as^ = 1/4 I 1 -2(g+) {S^)-{S^f- {sA . (39) 

cross ladder / 

The ladder terms oscillate quickly and are responsible for the fringes in Fig. 7. Similarly, 
the dispersion of Sy can be decomposed into the cross and ladder terms 



as, = 1/4 I 1 - 2 (5+) {S+f + {S^)\ I . (40) 



ladder 



The quickly oscillating ladder terms cancel out in the sum of the two dispersions 

^s. + ^5, = 1-{S^){S-). (41) 

The sum of the electronic dispersions up to a constant and a sign is nothing but the quasi- 
classical contribution 

£^quasici = hn{s+){s.) (42) 
17 



to the total electronic energy 

= hQ{S+S^). (43) 

Note that the corresponding sum of vibrational dispersions gives the quantum contribution 
to vibrational energy, Eq. (33). The difference between -Equasici -E^tai then is the purely 
quantum contribution to the electronic energy 

Substitution of the commutator [S^, Sy] = —iS^ and sum of electronic dispersions Eq. (41) 
into the general formula (28) gives the relation (37) between the electronic dispersions in an 
equivalent form 

crs^ + as, + (S,) > 0, (45) 

with Sz = S+S^ — 1/2 and Eq. (41). Eqs. (37) and (45) imply that the purely quantum part of 
the electronic energy is always greater than zero, -E'quant — 0- Over time the electronic energy 
flows between its quasiclassical and quantum parts, as shown in Fig. 7. The quasiclassical 
contribution to the electronic energy provides an envelope for the electronic dispersion. 

The electronic and vibrational dispersions show interesting similarities that are summa- 
rized in Table I. The fringes in the electronic dispersion. Fig. 7, appear at the same phase 
and frequency as the fringes in the vibrational dispersion. Fig. 5. The difference is in the 
sign of the deviation of the electronic and vibrational dispersions from their bounds, cf. the 
signs in inequalities (30) and (37). The fringes in the vibrational dispersion are directed 
upward, while the fringes in the electronic dispersion are directed downward. The electronic 
dispersion is bounded from above in contrast to the vibrational dispersion, which is bounded 
from below. 

Figure 8 displays the sums of dispersions for the vibrational and subsystems in comparison 
to the real valued electron-vibrational correlation 

^(A.) = {(/3)-(a)(5+)-(a+)(5+)}, (46) 

directly related to the auxiliary operator /3, Eq. (15). The correlation oscillates with the same 
low frequency as the electronic and vibrational dispersions. Rigorously, the combination of 
operators Eq. (46) can naturally appear as a sum of dispersions cta + cb for cross-operators, 
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composed of both fermionic and bosonic parts A = S+ + a, B = i (5*+ — a) and their 
conjugates. For the initial condition with no correlation between the fermionic and bosonic 
modes, and zero horizontal spin projections (5*+) I^^q = 0, the expectation value of Eq. (46) 
starts at zero and oscillates around zero. Up to zero point value and scaling constants, the 
dispersions of electronic subsystem and the correlation Eq. (46) behave as sine and cosine 
functions. 

When the population inversion stalls. Fig. 4, the electron-phonon system is in a state 
with large uncertainty. For both electronic and vibrational subsystems, states with larger 
uncertainties provide larger purely quantum energy contributions. As shown in Fig. 5, the 
dispersion associated with the vibrational subsystem increases, in accord with the uncer- 
tainty relationship (30). In contrast to the vibrational dispersion, the dispersion associated 
with the electronic subsystem decreases as seen in Fig. 7, in accordance with the uncertainty 
relationship (37). The electronic subsystem approaches the state of minimal uncertainty 
when inversion oscillation stops. 

The unusual behavior of the dispersion associated with the electronic subsystem can be 
understood with the concept of pure and mixed states. An isolated electronic system is 
characterized by either a normalized wave function {tp) = cq |0) +ci or a density matrix, 
or quantum mechanically averaged values of electronic operators. In an entangled electron- 
phonon system the electronic part is a statistical average over different realizations of the 
vibrational subsystem. The statistical mixture of states cannot be represented through a 
wave function, and is characterized by a reduced density matrix or both statistically and 
quantum mechanically averaged values of electronic operators. 

The density matrix of an isolated electronic system is determined by the products of the 
wave function coefficients pe\ = \ip) {ip\ = cjci |1) (1| + cJcq |1) (0| + CqCi |0) (1| + CqCq |0) (0|. 
The reduced density matrix of the electronic subsystem of the electron-phonon system can 
be specified by expectation values of electronic operators averaged over the state of the whole 
electron-phonon system pei = {S+S^) |1) (1| + {S+) |1) (0| + |0) (1| + {S^S+) |0) (0|. The 
trace of the squared density matrix of an isolated electronic system is a complete square and 
always equals unity Trp^j = (c^Ci -|- CqCq)^ = 1, due to normalization of the wave function. 
The trace of the reduced electronic density matrix squared, often referred to as fidelity, takes 



19 



the form 



Trp2 = 1 + 2(S+)(S_) - {S+S^){S^S+) < 1, (47) 

and is less or equal one. The last term is nonnegative, because it is a product of nonnegative 
populations {S±Szf) > 0. The second term is zero if the averaging over the phonon state 
gives no correlation between |0) and Trp^j = 1 means that the electronic subsystem 
exists in a pure state. Trp^j < 1 implies that the electronic subsystem is in a mixed state. 
The distinction between pure and mixed states is best illustrated with the Bloch vector 

R^{{S.),{S,),{S.)). (48) 

The square of the Bloch vector is related to Trp^j by a linear transformation given in 
Table II. The Bloch vector (48) connects the origin with the point in a three-dimensional 
space specified by the expectation values of the projections of the spin-1/2 operator onto x, 
y, z axes. 

= {S^f + {Syf + {S^ < 1/4. (49) 

The square of the Bloch vector can be decomposed R^ = {S^)"^ + R^^ into the vertical {S^) 
and horizontal components 

R% = {S^f + {S,f = (S^) . (50) 

The X, y, and z components of the Bloch vector of an isolated electronic 

system are {cIcq + CqCi) /2, (c^cq - CqCi) /2i, (c^ci - CqCq) /2 yielding R"^ = 

i [(cjci)^ + (cgCo)^ + 2ciCiCgCo] = 1/4. The squared length of the Bloch vector of 
the electronic subsystem of a coupled electron-phonon system is 

R^ = 1/4+ (5+) - {S+S^){S^S+) < 1/4. (51) 

Here, the last term is always nonnegative, and the second term is zero if the correlation 
between and 1 is dephased by averaging over the vibrational states. R^ decreases to 
zero when the electronic subsystem is equally mixed between the two states, {S+S-) = 
(S^S^) = 1/2, the last term compensating the 1/4 constant. The relationships between the 
square of the Bloch vector and the trace of the squared density matrix are summarized in 
Table II. As shown in Fig. 8(c) for electronic subsystems with equal populations of both 
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electronic states {Sz) = 0, {S+S-) = (S'_5'+) = 1/2, the change of the squared Bloch vector 
from to 1/4 corresponds to the change of fidehty from 1/2 to 1. Fidehty expressed in 
terms of the electronic dispersions shows that mixed electronic states correspond to larger 
values of the electronic dispersions. Pure electronic states yield minimal uncertainties of 
electronic operators, minimal values for the electronic dispersions and the maximal length 
of the Bloch vector R=l/2. Alternatively, mixed electronic states yield larger electronic 
operator uncertainties, larger electronic dispersions and shorter Bloch vectors. 

The end of the Bloch vector glides from the south pole of the Bloch sphere to the north 
pole then back to the south pole and so on. During this revolutions, shown in the lower 
panel of Fig. 9, the length of the Bloch vector decreases, as the electronic subsystem becomes 
a statistical mixture. The mixing is maximized at ut ~ 1.5 — 2.0. For times ut > 2 the 
Bloch vector becomes larger and the electronic subsystem returns closer to a pure state. The 
length of the Bloch vector never recovers its initial value during the studied time period. In 
the middle of the silent period the Bloch vector is almost at maximal length ~ 1/2. Its 
z-axis projection is zero (S^) = in accordance with Fig. 4. Thus, during the silent period 
the Bloch vector revolves equatorially (nutates) around z-axis. 

The analysis of fidelity shows that the electronic subsystem successively passes through 
the following sequence of states. Initially in a pure state, the electronic subsystem quickly 
becomes a mixed state with a small fidelity. This occurs within one-eighth of the revival 
period, defined by f^revivai — g/y/lo^^- Then, the fidelity rises and the electronic subsystem 
again closely approaches a pure state. The fidelity is maximized at exactly one half of the 
revival period, t = n/Qj-cvnai- 

During the silent period when the electronic population stalls, the oscillation shifts to 
the vibrational and electronic dispersions and (S^), (Sy) expectation values that preserve 
the phase of the oscillation and are maximized exactly in the middle of the silent period. 
Generally, the electronic and vibrational dispersions and electronic coherences {S^), {Sy) are 
maximized at half integer revival times trcvivai = 27r/f2revivai = Svr^/To/^f, while the electronic 
inversion (Sz) and vibrational classical-like (p) and (q) are maximized at integer revival 
times. 



21 



V. TIME SCALES 



The coupled electron-phonon dynamics in the molecular dimer forms a hierarchy of time 
scales. The electronic and vibrational coherences (^x) and (g), as well as electronic pop- 
ulation {S2) and energy (a+a) oscillate with different frequencies that vary by order of 
magnitude or more. Each frequency component can be associated with a distinct physical 
process. The following is the hierarchy of the time scales in increasing order. The quickest 
time scale tj corresponds to the oscillations of the expectation values of the vibrational co- 
ordinate (g) and electronic coherence (^x), (•S'y). The second time scale tjj characterizes the 
inversion of the electronic population {3^)- The population inversion time depends on the 
initial displacement of the vibrational wavepacket. Larger displacements yield quicker pop- 
ulation transfer. Relaxation and dephasing of the population inversion occurs on a slower 
time scale tm. Over this time interval the electronic state approaches quasi-equilibrium 
with equal partitioning of the population between the two states. The slowest component 
tjv is associated with the silent period and determines the revival of the oscillation in the 
electronic state populations. The tjv time also determines the time scale of the relaxation 
of the expectation value of the vibrational coordinate (q) to its equilibrium position. The 
electron- vibrational system enters and leaves the state with the maximal uncertainty on the 
tjv time scale. 

The time scales of the faster processes / through III can be estimated from the ap- 
proximate analytic expressions (20)-(21) derived in this work. The slowest time scale IV is 
estimated by Eberly et al.^^ and is reproduced in our numerical simulations. The slowest 
time scale IV describing the silent period and revivals of the oscillation in the electronic 
populations does not appear as a separate time scale in the analytic expressions (20)-(21) 
and coincides with the time scale /// of the relaxation and dephasing of the oscillation. The 
revivals in the electron-phonon system can take longer than the dephasing time and produce 
a silent period. According to Eberly et al.^^, tjv is determined by the difference in the Rabi 
flopping frequencies of the two neighboring and most populated levels of the oscillator. The 
following are the rates associated with the four time scales 



t 
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Even for this simple two dimensional system isolated from environment, the time taking the 
system to return to its initial state, the Poincare recurrence time can be very long. The 
quasi- irreversibility feature appears due to discreteness of the oscillator states and determines 
the character of the long-time dynamics. After each revival period the state of the system 
is farther away from the initial state. For most real systems with large number of degrees 
of freedom even the first revival is barely accessible. According to Figures 2 and 5, after 
a long period of time the expectation values of the population inversion and vibrational 
coordinate lose the revivals and approach quasi-equilibrium. Figure 7 indicates that some 
electronic coherence may persist for a long time. Based on the above analysis it is reasonable 
to expect that after a sufficiently long time interval the system will be found in a state of 
large uncertainty with the first order dynamical averages performing small amplitude chaotic 
motions about the equilibrium values. 

VI. DISCUSSION 

Few simple quantum models allow for analytic treatment: two interacting harmonic oscil- 
lators as well as two interacting two-level systems are exactly solvable. A harmonic oscillator 
interacting with a two-level system cannot be solved exactly by use of a finite number of 
variables. Known are the infinite series solutions of the eigenstate method and finite approx- 
imations by the QHD method, presented here. The infinite set of wave function coefficients 
in the equations of motion in the eigenstate solution is equivalent to the infinite chain of 
QHD equations. 

The dynamics of a massive particle in a harmonic potential, undergoes the same evolution 
as a quantized mode of a massless electromagnetic wave. The relaxation of the population 
inversion is irreversible in general, so that the dynamics can be referred to as reversible 
decoherence. An interesting issue is observed in Fig. 6(a) at ujt/{2'K) = 9.5. The dispersion 
of coordinate closely approaches its initial value. The states of the electronic and vibra- 
tional subsystems become correlated with time. This nontrivial behavior of the averages is 
associated with formation of entanglement. 

Applications and extensions of this work cover three directions: theoretical, chemical. 
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and spectroscopic. From the theoretical point of view the results can be improved by using 
higher order closure relations. The necessity to take into account zero-point energy effects 
and get a better description of the silent period suggests an extension of the calculation 
scheme by including equations of motion for new variables, in particular, the dispersions of 
products of nuclear and electronic operators. Variables which follow state-specific positions, 
electronic phase, and vibrational phase acquired by the dimer with time can be included 
into the calculation scheme for application to echo experiments. For instance, the state- 
specific positions are (g {S^ ± 1/2)). The relative electronic phase can be defined with {S^), 
(Sy). The vibrational phase is specified by a variable, analogous to the 7 term in the Heller 
work on Gaussian wavepackets^^. Some of these new variables and equations of motion 
appear naturally while describing the coupled electron-phonon dynamics with the Fokker- 
Planck equation of motion for the normally ordered multidimensional generating function, 
comprising both electronic and nuclear degrees of freedom^^. 

An important chemical application of the presented method of calculating the coupled 
electron-phonon dynamics is the description of energy and charge transfer in large molecular 
aggregates playing important biological and technological roles, e.g. in the light harvesting 
complex and molecular solar cells. The dynamics of electronic excitations is coupled in 
large systems to many inter- and intramolecular nuclear degrees of freedom as well as to 
vibrations of environment, whose dynamics are extremely difficult to follow on the quantum 
level. The QHD method used here follows the quantum properties using few expectation 
values to characterize each degree of freedom, and describes the exciton and electron transfer 
dynamics in a simple way. The QHD method is directly applicable to the quantitative 
description of transfer process in large chemical systems. 

The method of this paper is also able to describe 3 pulse-, 4 pulse-, and correlated 
photon echo spectroscopic experiments with molecular systems by addition of new variables. 
In femtosecond echo experiments with ensembles of disordered dimers^'^^, the evolution of 
phase for each dimer in the ensemble determines the re-phasing of the ensemble net transition 
dipole and, consequently, the photon echo signal. In particular, for the three-pulse stimulated 
photon echo technique, the phase acquired during the population period between pulses II 
and HI determines the time delay for appearance of the stimulated echo. The echo signal 
determined by the relative phase between the ground and one of the excited electronic 
states, carries information about the population transfer between the optically active excited 
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states. It is expected that oscillations in the third-order echo peak shift, modulated by 
the damped evolution of the excited states population {Sz), will be reproduced by the 
current method. In a four pulse femtosecond experiment, the fourth pulse serves for the 
heterodyning: The excited state wavepacket is projected onto the ground state wavepacket. 
The created interference population in the excited state affects the net fluorescence from this 
state. The time of the signal is determined by the acquired electronic phase and position of 
the exciton transfer wavepacket described by the {q{Sz ± 1/2)) variable. 

A novel correlated spectroscopic technique, based on a simultaneous irradiation of a sam- 
ple by pulses of optical (electronic) and IR (vibrational) frequencies can yield new informa- 
tion about the dynamics of the electron- vibrational correlation c.f. (a) and {(3), calculated 
in this paper. We propose to excite and probe by correlated infrared and optical pulses. 
The scenario is analogous to the 2D IR spectroscopy used to get the coupling between the 
vibrational modes with close frequencies^" . For the case of symmetric modes, the dynamics 
of coupling between optically active low frequency vibrational and high frequency electronic 
modes can be revealed by 2D Raman. For non-symmetric modes direct and independent 
excitation of the vibrational degree of freedom with an IR-pulse and the electronic degree 
of freedom with an optical pulse is required. The technical difficulty of this experimental 
method is establishing the phase relations between the IR and optical excitations. The ex- 
perimental data measured by such technique can be modeled by the method developed in 
this paper. 

VII. CONCLUSIONS 

Quantum dynamics of an electron-nuclear excitation in a model dimer has been analyzed 
in this work in the limit of small nuclear reorganization. The relaxation behavior in the model 
has been obtained by a reduced description that, by simplicity, is applicable to large systems. 
The approximate approach quantitatively reproduces the correct behavior of the model. It 
is also demonstrated that even in the absence of the heat bath the system approaches a 
mixed state corresponding to a quasi-thermal equilibrium, by destructive interference of 
wave function coefficients. 

The influence of the vibrational mode onto the electronic population dynamics is ana- 
lyzed by solving a chain of coupled differential equations for population inversion, electron- 
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vibrational correlation, etc. The derived analytic result describes oscillation and relaxation 
of electronic population at short times. The approximate solution correctly represents the 
first several population inversions and the decay in the amplitude of the oscillation of the 
electronic population. The approximation works better for the quasiclassical regime with 
large vibrational amplitude and weak coupling. 

Detuning of the electronic and vibrational frequencies slows down the relaxation of the 
electronic population oscillation. An increase of the initial vibrational energy qualitatively 
changes the evolution of the electronic population from continuous oscillation to beats that 
reflect relaxation. The rate of the relaxation drops with increasing vibrational amplitude. 
It is found that several dynamical variables oscillate with the same rate and relax at the 
same time. Quantum expectation values, such as the vibrational coordinate dispersion 
not included into the primary calculation scheme, preserve the dynamics during the silent 
period, when a single vibrational mode acts as a quasi heat bath and induces electronic 
population dephasing over a finite time. The relaxation of the electronic population and 
a noticeable growth of the vibrational coordinate dispersion originate from splitting of the 
vibrational wave function into two wavepackets that oscillate independently and correlate 
with different electronic states. The original Gaussian wavepacket transforms into a non- 
classical state with large width. The vibrational energy flows between its quasiclassical and 
purely quantum components. Relaxation of the oscillation of the overall vibrational energy 
is compensated by growth in the oscillation of the vibrational coordinate dispersion. 

The vibrational energy contains more quasiclassical energy at the initial and revival times 
and less quasiclassical energy during the silent time interval between the revivals. At the 
time of revivals the vibrational subsystem comes into a quasiclassical state. By analogy, 
the electronic energy flows between its quasiclassical and purely quantum components. The 
electronic energy contains less quasiclassical contribution at revival times and more quasi- 
classical contribution during the silent period. The interaction with the vibrational mode 
drives the electronic subsystem into a mixed state and reduces the length of the Bloch vector. 
In the middle of the silent interval, the electronic subsystem comes into a coherent state. 
This qualitative process is reflected in the time evolution of the squared length of the Bloch 
vector, the sum of the electronic dispersions, and the fidelity of the electronic subsystem, 
which are all linearly dependent. 

Time dependence of the sums of the electronic and vibrational dispersions, and the 



26 



electron-vibrational correlation oscillate with a low frequency corresponding to the rate 
of relaxation and dephasing of the fast first order dynamical variables. The extrema of the 
sums of the electronic and vibrational dispersions correspond to the fastest change in the 
electron-vibrational correlation and vice versa. Oscillations of the dispersions of the elec- 
tronic and vibrational operators occur with the same frequencies. The correlation between 
the electronic and vibrational subsystems grows, oscillates, relaxes, and reappears similarly 
to the population inversion. Work in progress extends the reported results toward more 
accurate calculations and descriptions of molecular charge and excitation transfer observed 
in ultrafast optical experiments. 
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TABLE I: Electronic and vibrational energies partitioned into quantum and classical contributions. 
The energies are given in the ladder, phase space and dispersion operator representations. In both 
electronic and vibrational subsystems the energy flows between the quantum and quasi-classical 
contributions. 



Energy 



total 



quasiclassical 



pure quantum 



Vibrational 



ladder operators a~^, a 
phase space operators p, q 
dispersions 
equation 



(a+a) + 1/2 {a+){a) 



2m 



(34) 



2m ^ 2 



(35) 



(a+a) - (a+)(a) 



2m 



+ 



1 „ I mcj^ „ 
2m '^9 2 

(30), (33) 



Electronic 



ladder operators, 5"+, {S^S-) (5+)(S'_) (5+5_) — (S'+)(S'_) 

projection operators, 5^, 5y, 5, (5,) + 1/2 {S^)^ + (5y)2 (5,) + 1/2 - ((Sx)^ + {Sy)^) 



dispersions 
equation 



{crs^ + CFSy) (5'z) + (crSx + CfSy) 



(43) 



(41) 



(37) 
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TABLE II: Key characteristics of the electronic subsystem: Bloch vector squared (i?^), pure quantum energy {E^^^^^), and fidehty 
(Trp^) are represented by the expectation values of the ladder, projection and dispersion operators. These characteristics obey the 
equalities for pure states and the inequalities for mixed states. 

i?2 < 1/4 > Trp^<l 

ladder, S+, 5_ 1/4 - (5+5_)(5_5+) + {S+S-) - 1 + 2(5+)(S_) - 2{S+S-){S-S+) 

projections, S^, Sy, (S^)^ + (Sy)^ + {S,f (S,) + 1/2 - (S^)^ - (Sy)^ 1/2 - 2{S,f + 2 {{S^)^ + (Syf) 

dispersions 1/2 + {S,f - as^ - as, {S,) + as^+as, 1/2 - 2{S,f + 1 - 2 {as^ + as,) 



Figure Captions 

FIG. 1 Electronic states and state-specific potential energy curves as functions of the 
vibrational coordinate of the dimer described by the Hamiltonian (12), appearing in the 
theory of molecular excitons and also known as the Marcus model. The electronic states 
are separated by the energy gap e = Ml. The equilibrium positions in the vibrational 
coordinates of the two electronic states are displaced by the value specified in Eq. (11) and 
related to the nuclear reorganization (a + a*)/2. 

FIG. 2 The expectation value of the population inversion operator Sz as a function of 
time computed exactly, solid line, and approximately by Eq. (20), dotted line, (a) illustrates 
the quasi-classical case of a weak electron-phonon coupling g = 0.025 and a large amplitude 
vibrational motion {a'^a){t = 0) = 49. (b) corresponds to a deep quantum case of a strong 
coupling g = 0.25 and a smaller amplitude vibration {a^a)(t = 0) = 9. (5'z)^^q = —1/2 
and (5 = in both cases. Good agreement between the exact and approximate solutions is 
observed during the first several periods of the population inversion. For both types of initial 
conditions, the approximate solution correctly reproduces the dephasing envelope describing 
the decay in the amplitude of the oscillation of the electronic populations. 

FIG. 3 Population inversion and relaxation tZ^ rates, Eq. (23), labeled by "+" 
and "— ", as functions of (a) detuning 6, and (b) amplitude a, Eq. (11). The remaining 
parameters are kept constant at quasi-classical {a~^a) = 49, g = 0.025 (solid line) and 
quantum (a+a) = 9, g = 0.25 (dashes) values used in Fig. 2. The dependence of the inversion 
relaxation rate on the amplitude of vibrational motion displays branching in (b). The 
inversion t'^'^ and relaxation tZ^ rates coincide for small values of the vibrational amplitudes 
a. The branches split at a = 1. 

FIG. 4 Evolution of the sum of the population inversion 2(5*2) and the electron-phonon 
correlation The purely real 2{Sz) and imaginary a^^{(3) functions shown in (a) are 

dynamic conjugates in the sense of Eq. (16) with n = 0. The absolute value \ {2Sz + a~^l3)\ 
is shown in (b). The initial conditions and parameters are {ci~^cL)^^Q = 9, (5'z)^^q = —1/2, 
g = 0.25, 6 = 0, corresponding to Fig. 2(b). The population inversion S^, and the 
electron-phonon correlation f] oscillate and dephase with the same rate. The amplitude 
I (25*2 + d~^P) \ displays smooth dephasing dynamics without the quick oscillations. 
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FIG. 5 Key characteristics of the vibrational subsystem, (a): Dispersions Ug, ap of 
the vibrational coordinate q and momentum p are shown by the solid and dotted lines, 
Eqs. (31) and (32), respectively. The sum + Eq. (33), shown by the dashed line starts 
from the minimum uncertainty width and smoothly increases to its maximum value that 
is equal to the total energy of the system. The sum of dispersions serves as an envelope 
for the individual dispersions. The slow increase of the dispersion is due to the splitting 
of the original Gaussian wavepacket into two Gaussians that oscillate independently and 
pass through each other leading to the rapid oscillations of o"g and ap. (b): Shown are the 
total energy of the vibrational mode Eq. (34), solid line, together with its quasi-classical 
contribution, Eq. (35), dots, and the purely quantum contribution, Eq. (33), dashes. The 
vibrational energy flows between its quasiclassical and purely quantum contributions. At 
t = the quantum contribution equals zero-point energy, and the rest of the vibrational 
energy is quasiclassical. When the total vibrational energy stops oscillating, it is purely 
quantum. When the oscillation resumes, the vibrational energy is dominated by the classical 
contribution. The sum of dispersions in (a) is an alternative representation of the quantum 
energy in (b), dashed lines. The initial conditions and model parameters are the same as 
in Fig. 2(b). Note, that the phase of the vibrational energy oscillation in this flgure is 
opposite to the phase of the electronic population in Fig. 2(b), due to conservation of the 
total electron-vibrational energy. 

FIG. 6 (a) Time slices of the vibrational wavepacket in coordinate representation, (b) 
the expectation value of the vibrational coordinate q, (c) and the expectation value of 
the population inversion (S^) + 1/2. The initial conditions and model parameters are taken 
from Ref. [59]. Transfer of population between the two states affects the wavepacket dynamics. 
The original Gaussian wavepacket splits into two wavepackets correlated with each electronic 
state. As a result of the splitting the oscillation of the average vibrational coordinate decays. 
Most interestingly, the oscillation of the electronic population decays as well. The decay in 
the oscillation of the electronic population, also known as dephasing, typically occurs by 
coupling to a heat bath. Remarkably, a single vibrational mode acts as a quasi-heat bath 
and induces population dephasing over a flnite time. 

FIG. 7 Key characteristics of the electronic subsystem, cf. Fig. 6 for the vibrational 
subsystem, (a): Dispersions as^, (JSy of the conjugate electronic operators and Sy are 
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shown by the sohd and dashed hnes, Eqs. (39) and (40), respectively. The sum as^ + o"Sy — 
1/4, Eq. (41) shown by the thick dot-dashed hne starts from the maximum value of 1/4, 
corresponding to a pure state, and varies more smoothly than the components as^, (JSyi 
creating an envelope. In contrast to the vibrational dispersions. Fig. 5, which increase due 
to the electron- vibrational interaction, the electronic dispersions decrease, (b): Shown are 
the total energy of the electronic subsystem {S+S-), Eq. (43), solid line, its quasiclassical 
contribution, Eq. (42), dots, and the purely quantum contribution, Eq. (44), dashes. The 
electronic and vibrational energies, this figure and Fig. 5, oscillate with opposite phases 
so that the total electron- vibrational energy remains constant. The electronic energy flows 
between its quasiclassical and purely quantum contributions. At t = both quantum and 
classical contributions equal zero since system is prepared in the state (5'z)(0) = —1/2 During 
the silent period, when the population inversion stops oscillating, the electronic energy 
contains the largest portion of the classical contribution, while the vibrational energy is 
purely quantum. The sum of dispersions in (a) is an alternative representation of the classical 
component of the electronic energy in (b), up to the linear transformation of Eq. (41). The 
initial conditions and model parameters are the same as in Fig. 2(b). 

FIG. 8 Slowly oscillating combinations of the expectation values. The initial conditions 
and model parameters are the same as in Fig. 2(b). (a). Time dependence of the sums 
of the electronic and vibrational dispersions, solid and dashed lines, respectively. Surpris- 
ingly, the sum of the electronic dispersions, Eq. (41), gives the quasiclassical contribution 
to the electronic energy, Eq. (42), while the sum of vibrational dispersions gives the pure 
quantum contribution to the vibrational energy, Eq. (33). (b): Imaginary part of the 
correlation between the electronic and vibrational subsystems. The expectation value of 
the correlator Eq. (46) is purely imaginary, its magnitude is displayed on the y-axis. The 
electron-vibrational correlator starts at and evolves on the slow time-scale of formation 
and dephasing of the electron-vibrational entanglement. The extrema of the sums of the 
electronic and vibrational dispersions in (a) correspond to the fastest change in the correla- 
tion in (b) and vice versa. The characteristics displaced in the figure are expressed in several 
ways in Table I. 

FIG. 9 Trajectories of the Bloch vector, Bloch vector squared, and a relationship be- 
tween the key characteristics of the electronic subsystem, (a): Bloch vector, Eq. (48), as a 
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function of time during the first one and a half population inversions, for the same initial 
conditions and model parameters as in Fig. 2(b). The Bloch vector is defined by the three 
spin projections of the electronic two-level system (S^), {Sy), {S^) and is shown together 
with a Bloch sphere of radius 1/2. At t = the Bloch vector connects the origin and the 
south pole of the Bloch sphere. The Bloch vector revolves from the south pole to the north 
pole and changes in length, (b): The squared length of the Bloch vector as a function of 
time. -R^ = 1/4 at t = 0, decreases at < u;t < 2, and rises to a local maximum at ut ^ 6. 
(c): Key characteristics of the electronic subsystem as functions of the horizontal projection 
of the Bloch vector, Eq. (50), for {S,) = 0, {S+S^) = {S^S+) = 1/2. Shown are the Bloch 
vector length squared Eq. (51), solid line, the sum of dispersions for the electronic subsystem 
Eq. (41), dot-dashes, the quantum contribution to the electronic energy Eq. (44), dots, and 
fidelity Eq. (47), dashes. The right edge, (^x)^ + {Sy)'^ = 0.25, corresponds to pure states 
characterized by maximal fidelity, Trp^j = 1 and = 1/4. The left edge, + (•S'y)^ = 0, 
describes the mixed state with minimal fidelity, Trp^j = 1/2 and = 0. Interaction with 
the bosonic mode drives the electronic mode from a pure to a mixed state. 
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